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THE VANISHING VISCOSITY LIMIT IN THE PRESENCE 
OF A POROUS MEDIUM 


CHRISTOPHE LACAVE AND ANNA MAZZUCATO 

Abstract. We consider the flow of a viscous, incompressible, Newto¬ 
nian fluid in a perforated domain in the plane. The domain is the ex¬ 
terior of a regular lattice of rigid particles. We study the simultaneous 
limit of vanishing particle size and distance, and of vanishing viscos¬ 
ity. Under suitable conditions on the particle size, particle distance, 
and viscosity, we prove that solutions of the Navier-Stokes system in 
the perforated domain converges to solutions of the Euler system, mod¬ 
eling inviscid, incompressible flow, in the full plane. That is, the flow 
is not disturbed by the porous medium and becomes inviscid in the 
limit. Convergence is obtained in the energy norm with explicit rates of 
convergence. 


1. Introduction 

This article concerns the vanishing viscosity limit for incompressible, New¬ 
tonian fluids in a perforated planar domain when the size and distance be¬ 
tween obstacles vanish. By a perforated domain, we mean the exterior of 
a (finite) regular lattice of rigid particles. This problem can be viewed as 
simultaneously taking the limit of vanishing viscosity for flows in exterior 
domains and study the permeability of a porous medium under this flow. 

From a physical point of view, it is important to study fluid flow through 
a porous medium in the regime of small viscosity, which model for example 
flow of underground water. Away from boundaries, the effect of viscosity is 
negligible generically for such flows, but near walls, viscous boundary layers 
appear where large stresses and vorticity can develop. The layer can in turn 
destabilize, which lead to boundary layer separation and the formation of 
a turbulent wake. (We refer to |32j for an introduction to the theory of 
boundary layers.) 
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We say that the vanishing viscosity limit holds on the time interval [0,T] 
if the solutions of the Navier-Stokes equations converge to the solution of the 
Euler equations strongly in the energy norm, i. e., strongly in L“([0, T]; L^). 
It is not known whether the vanishing viscosity limit holds in domains with 
boundaries, not even in two space dimensions, if no-slip boundary condi¬ 
tions are imposed on the velocity. The Reynolds number is a dimensionless 
quantity, which is inversely proportional to the viscosity. When the charac¬ 
teristic local length scale near the boundary becomes small at a faster rate 
than the viscosity, one can expect that the local Reynolds number stays of 
order one, preventing the formation of a strong boundary layer, as was ob¬ 
served in m, where the vanishing viscosity limit was established in the case 
of one shrinking obstacle. Homogenizing the Navier-Stokes equations in the 
perforated (periodic) domain with Reynolds number below a critical value 
yields Darcy’s law as in the case of the Stoke’s system, while homogenizing 
the Euler equations gives non-linear filtration laws that depend on the rela¬ 
tion between the particle size and the characteristic velocity of the problem. 
In the intermediate regime of high Reynolds numbers, asymptotic analysis 
leads to consider the Prandtl boundary layer equations, which have recently 
been shown to be ill-posed mm, for the cell problem. (We refer to |31j 
and references therein, in particular [30], for a more detailed discussion. See 
also |3| for a recent result on stochastic homogenization.) In the case of 
the Darcy-Brinkman system, the equations for the boundary corrector are 
linear and passage to the zero-viscosity limit is possible HHiEg. 

It is therefore important to identify regimes where it is possible to neglect 
the effects of both the viscous boundary layers as well as the porosity of the 
medium and use the Euler equations in the full plane. In this article, we 
establish the convergence of solutions to the two-dimensional Navier-Stokes 
equations in the perforated domains to solutions of the Euler equations in 
the full plane in the energy norm, when viscosity, size of the particle, and 
particle distance vanish in the appropriate regime: namely, particle size 
must be smaller than particle distance, and a certain ratio of particle size to 
particle distance must be bounded above by viscosity. We consider the most 
difficult and interesting case of no-slip boundary conditions for the Navier- 
Stokes problem. We confine ourselves to flows in the plane for technical 
reasons. However, our methods could be adapted to the more challenging 
study of 3D flows, at least for the case of a fixed number of obstacles (see 
the related discussion in m)- 

The study of 2D incompressible flows in the exterior of a single shrinking 
obstacle, both for viscous and inviscid fluids, was carried out in the series of 
papers da na da He] (see also [7] for results in the periodic setting), where 
the case of obstacles diametrically shrinking to a point was considered. This 
is the geometric set up in our paper as well. 

Even in the case of flows in the exterior of a finite number of shrinking 
obstacles at fixed positions, extending the results in [Tl] is not straightfor¬ 
ward, as one needs to choose an approximation to the Euler solution which 
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is divergence free and satisfy the no-slip condition at the boundary. For one 
obstacle, it is sufficient to truncate the stream function associated to the 
Euler velocity. For more than one obstacle, we gain a good control on the 
norm of the approximation by truncating the velocity and by constructing 
instead an appropriate corrector to restore the divergence-free condition. 
When the obstacles are kept at a fixed distance from each other, our result 
applies to arbitrary configurations of obstacles. 

When the problem is studied in a porous medium, additional difficulties 
arise. In fact, the permeability of the medium in the limit depends on the 
relation between particle distance and particle size. This relation changes 
drastically from the viscous to the inviscid case (cf. [U [2] for the viscous 
case, and igiHliaEi] for the inviscid case). Therefore, a careful analysis of 
the flow in the perforated domain is needed to pass to the limit, in particular 
with respect to the choice of initial data for the Navier-Stokes system. In 
the regime considered in this paper, we take the initial data as well as the 
Euler correctors to be independent of viscosity. We will expand on this point 
further in Subsection o 

We close with an outline of the paper. In the rest of the introduction 
we describe in detail the perforated domain and recall the fluid equations: 
the Navier-Stokes equations modeling viscous flows, and the Euler equations 
modeling inviscid flows. 

We will compare the Navier-Stokes equations in the exterior of the obsta¬ 
cles, extended by zero, to the Euler solution in the whole plane. To this end, 
we will construct suitable approximations of the Euler solution in the exte¬ 
rior domain that satisfy the no-slip condition at the boundary and are still 
divergence free. In the case of one obstacle, the authors in |16j truncate the 
stream function, as both the stream function and the pressure are defined up 
to an arbitrary constant, and can then be chosen to be small on the obstacle. 
Smallness of the stream function and the pressure avoid a sub-optimal esti¬ 
mate on the gradient of the truncation. This construction does not extend 
to more than one obstacle. In Section we present a different construction, 
which utilizes a direct cut-off function on the velocity and a correction to 
restore the divergence-free condition. This construction was introduced in 
[22] and is based on Bogovskii operator. Section also contains the elliptic 
estimates that are needed to establish our main convergence result. We also 
recall some standard estimates for the Euler solutions. 

In Section]^ we discuss the convergence of the Navier-Stokes initial data 
to the Euler initial data under the conditions on e (the particle size) and 
(the distance between particles) given in Theorem |l.l[ using the analysis in 

Ej- 


Finally, in Section]^ we prove convergence of the Navier-Stokes solutions 
to the Euler solution n® by correcting the Euler solution thanks to the 
results in Section]^ In fact, we prove a more general stability result— see 

analogous to the main result in |16j . The 


the statement of Theorem 4.1 
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convergence follows from this stability result and the convergence of the 
Navier-Stokes initial data to the Euler data, discussed in Section 

Throughout the paper V will denote the gradient of a function with re¬ 
spect to its argument, while C will denote a generic constant that may 
change from line to line. We will represent a point in by a: = (xi,X 2 ) 
and denote the ball in (disk) with center x and radius R by B{x, R). 

Acknowledgments. A. M. would like to acknowledge the hospitality and 
financial support of Universite Paris-Diderot (Paris 7), Universite Pierre et 
Marie Curie (Paris 6), and the financial support of Paris City Hall and the 
Fondation Sciences Mathematiques de Paris, to conduct part of this work. 

The authors are also grateful to the anonymous referee for his valuable 
comments on the first version of this article. 


1.1. The perforated domain. We begin by describing in detail the geo¬ 
metric setup. By a perforated domain we mean the complement of finitely- 
many inclusions arranged in a lattice. In this article, we consider the case 
where all the inclusions have the same shape and the lattice is regular: 


( 1 . 1 ) 

where 

( 1 . 2 ) 


JCf j := zfj + eJC 


/C is a connected, simply connected compact set of 

/C C [-1, l]^ die E and 0 ek ■ 


The regularity of the lattice is needed only when the distance between the 
inclusions vanishes. If a fixed number of shrinking obstacles is considered 
instead, these can be arranged in an arbitrary fashion. 

In (1.1), zfj should be chosen such that the inclusions are disjoints. We 


will assume that the inclusions are at least separated by a distance 2ds'. 


(1.3) 


^i,j + 2(f — 1)(£ T c^e)) £ + 2(j — l)(e -|- ds)) 
=ie,e) + 2(e -h 4)(* - 1, J - 1), 


see Figure We are interested in the case for which d^ goes to zero as 
e —>• 0, hence the case for which the number of inclusions goes to infinity. 

In the horizontal direction, we distribute the inclusions on the unit seg¬ 
ment, i.e. we consider 


f = 1,... ,nf 


with 


nf = := 


- l + 2de - 
-2(e ds)- 


where [x] denotes the integer part of x. In the vertical direction, we charac¬ 
terize the geometry by introducing a parameter 0 < /U < 1. The case = 0 
corresponds to the case in which we only have obstacles arranged on a line, 
whereas the case Ai = 1 corresponds to the case in which we have obstacles 






VANISHING VISCOSITY IN POROUS MEDIUM 


5 


arranged in a square lattice. This configuration is achieved by considering 
horizontal lines of obstacles, that is, 


j = l,...,n| with 


n| : = 


{N‘ 


£\ll 


The intermediate case 0 < fi < 1 corresponds to a rectangular configuration 


o o 

9 o 

24 


o 

o 

o 


L= 1 


o 

o 

o 


Figure 1. Distribution of inclusions. 


of inclusions, such that in the limit e —?■ 0 the density of particles becomes 
infinite, but the height of the rectangle vanishes. For simplicity, we carry 
out the analysis only in the case that all the inclusions have the same shape 
/C, but similar results can be established when obstacles of different shapes 
shrink homotetically to zero. When the number of inclusions becomes infi¬ 
nite, we need to assume that the possible shapes are finitely many, in order 
to have uniform constants in elliptic estimates (see Subsection 2.1). The 


main reason for studying different arrangements of inclusions as e goes to 
zero is that the limiting equations may change depending on the geometric 
configuration (this is indeed the case for the viscous homogenization of the 
perforated domain mmy 

Throughout the paper, the fluid domain will be the exterior of the obsta¬ 
cles: 


(1.4) 




\(uu 


/C: 




i=lj=l 


For notational convenience, we suppress the dependence of the fluid domain 
on /i. 

The main purpose of this paper is to study the simultaneous limit of van¬ 
ishing viscosity v and the limit of vanishing e for the solutions of the Navier- 
Stokes equations in the exterior domain For different values of pi, we 
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determine a relation between v, £, and de such that the Navier-Stokes solu¬ 
tion in the perforated domain converges to the solution of the Euler equations 
in the full plane. 

The asymptotic of fluids outside obstacles that shrink to points is a very 
active area of research. For inviscid flows, Iftimie, Lopes Filho, and Nussen- 
zveig Lopes have treated in [T3] the case of one shrinking obstacle in the plane 
(extended to the case of one shrinking obstacle in a non simply-connected 
bounded domain by Lopes Filho in [26]), whereas Lacave, Lopes-Filho, and 
Nussenzveig-Lopes have studied in |23j the case of an infinite number of 
shrinking obstacles. These works consider also the case where the circu¬ 
lation of the initial velocity around JCf j is non zero, but do not provide a 
control on the distance between the holes. In the case of zero circulation 
(see (1.9)), Bonnaillie-Noel, Lacave, and Masmoudi prove in [6] that, when 
dg = solutions of the Euler equations in the exterior domain defined 
in (1.4) converge to the Euler solution in the full plane if a < 1 (more pre¬ 
cisely, if a < 2 — fi). Hence, in the limit the ideal fluid is not perturbed 
by the perforated domain if d^ = e“ with a < 1. For larger a, Lacave and 
Masmoudi establish in [2Tj that the perforated domain becomes imperme¬ 
able (e.g. when /r = 1, an impermeable square appears if a > 1). In the 
periodic setting (related to the case /a = a = 1 in the context of this paper), 
an homogenized limit for a modified Euler system was obtained in [251 EH- 

Concerning viscous fluid with fixed viscosity n, Iftimie, Lopes-Eilho, and 
Nussenzveig-Lopes considered again the case of one shrinking inclusion in 
m (see also |13j in dimension three). There is also an extensive literature in 
the homogenization framework. Following the pioneering work of Cioranescu 
and Murat [^ for Laplace’s equation, Allaire studied the homogenization of 
the Navier-Stokes equations in a perforated domain under different regimes, 
given in terms of relative ratios of de to e, in PEj. When 11 = 1 and 
de = l/Vl lne|, or when /r = 0 and de = 1/| lne| for the critical distance 
between inclusions, the limit problem is described by a filtration law of 
Brinkman type. If the distance is larger than this critical value, one recovers 
the solution of the Navier-Stokes in the full plane without any influence of 
the porous medium, and if the distance is smaller, one obtains Darcy’s law 
in the limit. We refer to Pi and references therein for a more detailed 
discussion of this point, in particular for results on viscous flow through a 
sieve. We remark that the dependence of the critical distance on e needed to 
feel the presence of the porous medium in the limit is markedly different in 
the viscous and inviscid cases. For example, for distances de satisfying e“ <C 
dg <C 1/-^1 lne|, the limiting system describing the homogenized viscous 
flow is different than the Navier-Stokes system, whereas there is no memory 
of the obstacles in the limit for ideal flows. As a matter of fact, we are 
able to treat the situation in which the distance is smaller than the critical 
distance (namely, d^J^/\ln£\ —)■ 0) and the inclusions are distributed along 
a line, i.e., /i = 0 in our setting, albeit with a less than optimal rate of 
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convergence. For Navier-Stokes, this case was not addressed by Allaire in 
p] (see Remark 4.2). 

It is therefore natural to investigate the concurrent limit e, u —)• 0. There 
are also important physical motivations for studying this problem, as dis¬ 
cussed in the introduction. (We again refer to [21] and m for a more 
in-depth discussion of the physical implications of these and our results.) 


1.2. Fluids equations and initial data. We consider the flow of a vis¬ 
cous, incompressible, Newtonian fluid in 11^ with classical no-slip boundary 
conditions on := We hence let denote the solution 

of the Navier-Stokes equations with (kinematic) viscosity coefficient u in 
and initial velocity Uq: 


(1.5) 


’ut’^ - + (n^'= • = 0, 

divu'^’^ = 0, 

= 0 , 

_n^:£(0, •) = n5. 


on (0, -l-oo) X 
on [0, -|-oo) X 
on (0, -l-oo) X 
on 


We take initial data depending only on e and not directly on viscosity u for 
reasons that will be clear below. However, this set up can be easily relaxed. 
The initial data is assumed only to satisfy the no-penetration conditions 
forced by the impermeable boundaries. By abuse of notation, we will often 
refer to the velocity u'^’^ as the Navier-Stokes solution, without mentioning 
the pressure 

We will also consider inviscid flow in the whole plane and denote by 
the solutions of the Euler equations with initial data uq: 


( 1 . 6 ) 


uf + (u^ • V)u^ -b Vp® = 0, on (0, -boo) x 
< divu^ = 0, on [0, -boo) x M?, 

u^{0, •) = uo, on M^. 


Our main result is a convergence result of to under suitable con¬ 
ditions on the relative strengths of u, £, and ds, assuming convergence of the 
respective initial data. Such a convergence can be achieved if Uq is chosen to 
be an appropriate truncation of the Euler initial velocity uq that is tangent 
to the boundary of and divergence free. Uq will be close to uq owing 
to the fact that the size of the obstacles is small. A standard approach to 
constructing the truncated velocity is to give the initial Euler velocity uq 
in term of a fixed initial vorticity independent of the domain (as e.g. in 
ini in o HSl ESI 123 Ea [Ml) ■ More precisely, let wq be a smooth function 
compactly supported outside the porous medium (i.e. compactly supported 
in \ ([0,1] X {0}) if p € [0,1) and in \ ([0,1]^) if p = 1), and define uq 
by 


1 


{x - y)-^ 
\x — yp 


wo(y) dy, 


(1.7) 


uo(x) = A:r 2 [wo](a;) : 
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where K ^2 denotes the Biot-Savart operator in the full plane. That is, uq is 
the unique solution of the problem 

(1.8) divMo = OinM^, curluo = wq in lim |tto(x)| = 0. 

X^OO 


It is well known (see e.g. m) that, with this initial data, the Euler sys¬ 
tem (1.6) has a unique solution . The properties of such a solution will 
be recalled in Subsection 


Remark 1.1. The condition that ujq be supported outside of the obstacles is 
used mainly in proving convergence of the initial velocities in L?' in Section 
which is then utilized in the energy estimates of Section In fact, to 
show that — uq belongs to one needs 



and the condition on the support of ojq guarantees that this condition holds 
uniformly in e. By comparison, in the different context of flow in a half¬ 
plane, Maekawa recently proved that the vanishing vioscosity limit holds, 
if the initial vortieity is (compactly) supported away from the boundary, for a 
time that is at least comparable to the distance of the support of vortieity from 
the boundary. Note, however, that in Theorem [m the vanishing viscosity 
limit is established for arbitrarily long intervals of time over which we cannot 
exclude that the support of the Euler vortieity intersects the obstacles. 


For any e > 0, there exists a unique vector field Uq defined on verifying: 
divug = 0 in curlug = wq in Uq • n = 0 on dO.^, 


(1.9) 




Uq ■ T ds = 0 for all i,j, lim |uo(^)l = Oj 


see m- The convergence of Uq to uq will be established in Section using 
the methods of [6]. As we will show in (3.1), the vector field Uq constructed 
in (1.9) belongs to L°°nL^’°°(n^). Hence, we can apply the result of Kozono 


and Yamazaki |21l Theo. 4] to conclude that the problem (1.5) has a unique, 
global-in-time, strong solution. 

We close by observing that neither the Euler nor the Navier-Stokes solu¬ 
tion have finite energy, nevertheless we will be able to show convergence, as 
n, e vanish, of u'^’^ — to zero in the energy norm. 


1.3. Main result. In this work, we are interested in studying the limit of 
vanishing viscosity —)• 0 at the same time as the limit e —)• 0, in which the 
obstacles shrink to points. In the case of a single obstacle, it was proved in 
[T6] that, for any given 0 < T < oo, there exists a constant Ci = Ci(/C, uq, T) 
such that, under the conditions £ < Civ, converges to in the energy 
space L°°((0, T), L^(M^)). In this result and in our main theorem, and 
Uq are extended by zero to the whole plane. Unfortunately, the arguments 
in m use in a crucial way the fact that we have only one obstacle: namely. 
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the authors utilize the stream function and choose a pressure for the 
Euler system that is zero at the center of /C, in order to construct a corrector 
for the Navier-Stokes solution in the fluid domain via truncation. Taking 
and small as the obstacle shrinks is needed to balance the growth 
of gradients of the cut-off function. Such a choice is not possible in the 
case of more than one obstacle, as there is only one degree of freedom for 
the Euler pressure and stream function in the whole plane. Inspired from a 
method developed in [22] , we will instead directly truncate the velocity field 
with good estimates on the truncation, which will allow to adapt the energy 
argument of m to establish the zero-viscosity limit. 

Due to the inherent difficulty with the vanishing viscosity in bounded do¬ 
mains, we are only treating the case where dg > e (without loss of generality 
we can assume that < 1), hence the total number of inclusions satisfies 

( 1 . 10 ) 

Our main result is the following theorem. 


Theorem 1.1. Given ujq G let the solution of the Euler equa¬ 

tions jg in the whole plane with initial condition uq (given in terms of uq 
in (1.7)j. For any e,u > 0, let d’^ > e and let he defined in (1.1)-(1.4). 
Let he the solution of the Navier-Stokes equations (1.5) in with ini¬ 
tial velocity Uq, given hy the unique solution of (1.9)J. Then, there exists a 
constant A depending only on JC such that if 

Au 




< 


ll‘^o||LinL°°( 

and if (jJo is supported in , then for any T > 0 we have 


(1-11) sup 11^2(0^) < Bt ’ 

o<t<T ^ 

where Bt is a constant depending only on T, ||a;o||2,iniui’°°(R2), and JC. 

Remark 1.2. i) In the case of a fixed number of obstacles shrinking 
homotetically to points (arbitrary located), the above theorem is a 
direct extension of the results in m, namely, under the condition 
e < Au (with A depending on the distance and the number of obsta¬ 
cles and on ||wo||LinL°“(]R2) j we obtain the convergence of the Navier- 
Stokes solution to the Euler solution with a rate of convergence of 
order y/v. 

ii) The fact that our Euler correctors (see Proposition \2.4^ are e-dependent, 
but not u-dependent, can he seen as a consequence of the fact that 
the local Reynolds number, built in terms of the characteristic size of 
the obstacles, stays of order one as v and e go to zero in the regime 
of Theorem o In this situation, the boundary layers are negligi¬ 
ble, as already observed in [T6| . In fact, the proof of Theorem M is 
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similar to the proof of Kato’s criterion, which implies the vanishing 
viscosity limit dZllMlEe]. 

iii) The additional factor ,fi+^uU 2 (|1.11D can he viewed as an effect of 

the homogenization of the porous medium. Using standard correc¬ 
tors as for the Laplace problem (see [Il[2l|8lE3]A one would obtain a 
smaller bound of the form (i +^)/2 / The limit of -y,/] lne| 

determines precisely when the so-called “strange term” in the homog¬ 
enization appears. However, this improvement is less relevant from 
an application standpoint, as discussed next, and would require a sig¬ 
nificant amount of additional technical work (cf. the discussion in 
Remark f.2). 

iv) Refined consequences of our main result arise by considering specific 

regimes for the ratio of particle distance to particle size, as discussed 
at the end of Subsection \n\ For instance, when p. = 1 (inclu¬ 
sions are uniformly distributed in a square) and d^ = If the 

critical distance for the homogenization problem, if the viscosity u 
satisfies Ce-sJ\ lne| < <C 1/| lne| — for example, for some 

a G (0,1) — our result gives convergence —>■ , whereas if 

viscosity is kept fixed as e —>■ 0 does not converge to the Navier- 
Stokes solutions in the full plane, as recalled already. Similar con¬ 
clusions can be drawn for smaller inter-particle distances, namely 
the theorem implies the convergence to n® if e d^, as in the case 
de = e“ with a < 1/3. When ^ = 0, the critical distance between in¬ 
clusions for the homogenization problem [2] is de = 1/| lne|, whereas 
in our result the limit holds when e <C dg' , as for the case dg = 
with a < 


2. Preliminaries 

We begin by presenting some basic elliptic estimates and estimates on the 
Euler solution that will be used throughout the article. 


2.1. Basic elliptic estimates. For the energy estimate in Section we 
will need to approximate a divergence-free vector field u defined in by 
divergence vector helds verifying the Dirichlet boundary condition on 
dLl^. To this end, we introduce an appropriate cut-off function as follows. 


As we consider the case where d^ > e, we deduce from (1.1)-(1.3) that 




U zF + e[-l,lY 


V(f, j) 


and that 

(2.1) [zF + e{-2,2f'^[^(^zl^^ + e{-2,2f'^ =% 'i{i,j)fi{p,q). 
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4){x) = 1, |x|oo < 3/2, 
^(x) = 0, |x|oo > 2, 


We let i/) > 0 be a smooth, cut-ofF function such that: 

( 2 . 2 ) 

where |x|oo = max, |xi|, and set 

(2.3) 


^p) = i-EE 

i=l j=l 


e ri^ 

^ ^ ^ X — 




X G 


Then, 0 < c/*" < 1 in i/)'" = 1 away from the obstacles and = 0 on a 
small neighborhood of the obstacles. More precisely, denoting 


(2.4) :=U(zT +£(-2,2)2 \e/c), 

hi 


we observe that Supp(l — 
nnmber of obstacles nfni, 
and conclude that 


(j)^) C Ai;. As we have the bound (1.10) on the 
we can easily estimate the Lebesgue measure of 


(2.5) 


eVp 

1 - +£||V(/>‘^||lp < Cp l<P<oo, 

de 


with Cp a constant independent of e. 

As (j)^u is not divergence free, the stream function ?/, satisfying u = 
can be truncated instead as in [HI m da uni, that is, one sets := 
V'*-((/>^i/;), which is divergence free and vanishes at the boundary. In this 
case, one needs to compensate the growth of the gradient of (ff as the obsta¬ 
cles shrink, which is of order 1/e. For a single inclusion shrinking to a point, 
which we identify with the origin, one can choose tp such that 1 /( 0 ) = 0 and 
then = (p^u -|- = 0(1). Of course, such a procedure cannot be 

applied if we consider more than one obstacle or if the obstacle is shrinking 
to a curve instead to a point. Here we present an alternate way to truncate 
divergence-free vector fields so that they have support in the complement 
of several disjoint obstacles. The following method was used in [52] to treat 
the case of an obstacle shrinking to a curve and it is related to the Bogovskh 
operator (see W)- 

For 1 < p < oo and for any e > 0 fixed, we define an equivalent norm 
II • 11^1,p in the Sobolev space by 


( 2 . 6 ) 


w, 


l,p .— 


sP 


Uaa + w^fWUA. 


i/p 


Lemma 2.1. Given 1 < p < oo, there exists a constant Cp depending only 
on p, such that the following holds: given any divergence-free vector field 
u{t, ■) G L°°(M^) for all t G [0,oo) and given any e > 0, the problem 

div h^ = • u, on [0, oo) x A^, 
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has a solution h^{t^ •) G Wq'^{A!.) satisfying 


\\h%t,-)\\^i.p < Cp\\V(j)^ ■u{t,-)\\LP(A,) Vt E [0,oo), 

Wit,-) - h%s,-)Wi.p < Cp\\V(j)^ • {u{t,-) -u{s,-))\\lp(^a,) yt,se [0,oo). 
Moreover, if E 17{Ae) then 


dh%t, •) 


dt 


w, 


1,P 


<Cp 


Vc/>^ 


du{t, •) 


dt 


LP{Ae) 


Later on in Proposition \2A\ we will employ this Lemma to construct a 
suitable corrector to the Euler solution that is supported away from the 
obstacles. Informally, extending h^ by zero, we readily verify that 7 := 
(j)^ u — h^ is divergence free, agrees with u away from the obstacles, and 
the 17 norm of 7 is of order ||rt||LP + e||V(/)^||Lp||rt||L°° (cutting the stream 
function would instead give a bound of order ||m||lp + ||V(/)^||lp||i/)||loo). 


Proof of Lemma 2.1 


(2.7) 


We start by introducing the bounded open set: 
U := (-2,2)2\/C. 


As U satishes the cone condition (see e.g. [9l Rem. III.3.4] for the dehnition). 
Theorem III.3.1 and Exercise III.3.6 (solvable thanks to Remark III.3.3) in 
[9] state that there exists Cp with the following property. Given any function 
f{t,x) such that fit,-) E LPfU), and Jjjfit,-) = 0 for all t, there exists a 
solution h{t, -) E Wq’^{U) of the problem: 

(2.8) div h = f, on [0, oo) x U, 


such that 


(2.9) II^(^e)IIw1'P([/) ^ C'p||/(t, •)||2,p([/) yte[ 0 ,oo), 

(2.10) \\h{t, -) - h{s, •)||ivi.p([/) < Cp\\f{t,-) - f{s, -)\\lp{u) Vt, s E [0, oo). 


Moreover, if ^ E 17{U) then 


dh{t, -) 
dt 


<Cp 

W^’P{U) 


dfjt, -) 
dt 


LP{U) 


Solutions are non unique, and the main difficulty addressed in |9j is how to 
obtain a solution h vanishing at the boundary and verifying the estimates 
with constant Cp depending only on U and p. The existence of such a 
solution follows from an explicit representation formula due to Bogovskh 

lU [5]. 

We next utilize this result to construct a solution h^ on the domain 

= U • 

hi 
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Let 11 be a given vector field on such that it(t, •) G for any 

t G [0, oo). For any i = 1,..., nf and j = 1,..., nl fixed, the function: 


fij{t, x) := e • u) {t, zfj + ex) 

is defined on the set U, is bounded (hence it belongs to U‘{U)), and satisfies 
for each t: 


[ = - [ div{(j)'^u){t,y)dy 

Ju ^ 




u ■ nds = - 


S Jr’S 


zf .+Ed(-2,2)‘2 £ Jzf ■+e{-2,2)^ 


(jj^u ■ n ds 
div It = 0, 


where we have used twice that u is divergence free on that (/>^ = 1 on 
zfj + £d{—2, 2)^, and that = 0 on zfj + edK. 

Therefore, we can apply Galdi’s results recalled above with / replaced 
by fij{t,x) to conclude that there exists a function hij on [0,cx)) x U such 
that hij{t,-) G Wq{U) is a solution of (2.8) with / = fij and satishes 
(2.9)-(2.10). Finally, we extend hij by zero on and we define: 

h^{t,x) := ( t. 




X — z: 




To finish the proof, we show that has all the properties stated in Lemma 
By definition, h^{t, •) G Wq'^{Ai,) for all t and we have 


2.1 


1 


div/i^(t,x) = -^(div/ijj) it 




X — z: 




1 


E/- 






X — z: 




1 




= V(j)^ ■ u{t, x), 


using that the obstacles are separated and SuppV(/>^ C + £^}- In 

the formula above, ,+^/j denotes the characteristic function of the set 
zfj + eU. Next, for any t, we compute 






2f,+£t/ 


-1^ 




, „ Izf .+£{/ 

2 ,J 


I 


-Vh 


t, 


X — z: 




dx 




t, 


X — z: 




dx 








<£^(7^^ / zfj + £x) dx 

<Cl\\Vct>^-u{t,-)rLHA,y 
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A similar calculation applies to derive the bounds on \\h^{t, •) — h^{s, ■ 
and \\dth%t,-)\\^i,p. 


□ 


Remark 2.1. In the proof of Lemma 2.1 above, we have used several times 
the fact that the obstacles are well separated (see Equation {2.1)). Therefore, 
our analysis can be applied only to the case d^ > Cs for some C independent 
of e. For distances d^ <C e, we would need to understand the behavior of Cp 
on domains of the form (—1,1)^ \ pK, as p —>■ 1“. 

In the same way, changing variables in the classical Poincare inequality, 
we obtain the following result. The proof can be found in |16t Lem. 3]. We 
sketch it here for the reader’s convenience. 

Lemma 2.2. Let x be a smooth cutoff function such that y = 1 in B{0,3) 
and X = 0 in il(0,4)”’. There exists a constant Ki depending only on K. such 
that for any u verifying 

XU G 

then 


( 2 . 11 ) 




where is given in Equation (2.4). 

Proof. We note that C il(0,1) for all e < 1. 

By the usual Poincare inequality, 

\\u\\L2(^^) < Ki\\Vu\\l2(^up 

for all u such that xu G Hq{B{0,4:) \ JC), where U is defined in (2.7). Then, 
from the definition of given in (|2.4|), it easily follow that 




< 


V / |n(rc)pdx = V / \u {z^ + ey)f e'^dy 
i,j ij 

J \eVu{zP+ey)\‘^ dy = s‘^Kl\\Vu\\l 2 ^^^y 


which ends the proof. 

We finish this subsection with a Sobolev-type embedding. 


□ 


Lemma 2.3. Let x be a cutoff function as in Lemma\2.2^ and let A^ again 


be given in (2.4). Then, there exists a constant K 2 depending only on K. 
such that for any u verifying xu G it holds 

II''^IIl4(a^) < ^/eK2\\Vu\\l2(^J^^y 

Proof. Bringing together the Sobolev embedding of H^{U) in L‘^{U) and the 
Poincare inequality, we have 

||lf||L4([/) < ^\\^\\m{U) < KK2\\Vu\\L2(^^y 

for all u such that xu G Rq(M^ \ fC). 
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Hence, a change of variables gives 

II«IIl4(a,) = (Z] / <'^( f \u{x)\'^dx) 

^^Jzf.+eU ^ ^ 


1/2 




- Z] (/ 1 '“ 1^ ' 

[ \sS/u{zfj + ey)\‘^ dy = eK^\\Vu\\l 2 (^A^y 
’J u 


□ 


2.2. Basic estimates on the Euler equations. We end these preliminar¬ 
ies by collecting some results on the Euler solution which will be used 
mainly in Section]^ Except for Proposition 2.4 the statements presented 


here are well known and complete proofs can be found, for example, in 
As mentioned in the Introduction, we give the initial velocity as uq = 
[wo], where K ^2 is the Biot-Savart kernel defined in(1.7) and the initial 
vorticity uiq £ C'“(M^). Hence, by a result of McGrath [29] there is a unique 
global strong solution of the Euler equations (1.6) in the full plane. 

The vorticity a; := curln^ verifies in a weak sense the transport equation: 

dtOJ + ■ Vuj = 0, 


which allows us to deduce the conservation of the norm of the vorticity: 

•)IIlp(]R 2 ) = I|wo||lp(E 2 ), Vt > 0, l<p<oo. 

By a classical estimate on the Biot-Savart kernel, we obtain that is 
uniformly bounded: 

(2.12) llri ||ioo^]]j+ xr2^ < ^II^IIl°°(r+xr 2)ll^ll£oo(-]]j+.2^i(']]j2p ^ C*ll^o||LinL°° 
and by the Calderon-Zygmund inequality, we infer that for all p G (l,oo): 

(2.13) ||Vn'®||ioo(]R+.2,p(R2p < C'p||a;||2,oo(R+.ip(R2p < . 

Hence, the nonlinear term ■ Vu® belongs to L°°(M+; n L^(M^)). Since 
p^ is a solution of Ap® = div(ri^ • Vn^), up to choosing p^ with zero mean 
value on H(0,2), it belongs to n 1E^’^(H(0, 2))) and we deduce 

from elliptic estimates and the Sobolev embedding that 

(2.14) p^ G L“(M+ X H(0, 2)). 

These estimates also imply that uf is in L“(M"*“; L^(H(0,2)). 

Since the Calderon-Zygmund inequality is not true for p = oo, there is 
no bound on ||Vtt®(t, •) ||x,oo(r 2 ) uniformly in time. However, the following 
well-known estimate |37| holds: 

(2.15) ||Vn^(t, •)IIl-(e 2) < Vt > 0, 

where Cq depends on ||a;o||vi/i,tx>- (See [20] for a discussion of results con¬ 
cerning the sharpness of the double exponential growth of ||Va;||i;,oo.) In 
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this paper, we do not investigate the minimal regularity of the initial data 
needed for our result to hold, as the vanishing viscosity limit is a singular 
problem even for smooth data. In particular, it is known that (2.15) holds 
under much weaker conditions (e.g. wq in the Besov space 135 ]). 

In the following proposition, we utilize Lemma 2.1 to construct a corrector 
to the truncated velocity (j)^u^ to meet the no-slip boundary conditions and 
the divergence-free condition, which will be employed for the energy estimate 
in Section m 

Proposition 2.4. Let coq G and let be the solution of the Euler 

equations (1.6) with initial condition uq, given in terms of ujq in (dZl);. For 
any e > 0, let and be defined as in (2.3)-(2.4). Then, there exists 
h£ g L“(]R+, such that, after extending h^ in by zero, and 

setting 

■= (jfu^ - h^^ 

we have that is divergence free and identically zero on . Moreover, 
there exists a constant Ke, dependent on but independent of e, such that 
for all £ > 0; 


(1) 

(2) imwv 

(3) \\Vh^\\E 

( 4 ) 


E ell 
U — U^Wlo 


41^)74 

]R+;L2(]R2)) < Ke jfdy 4 ,; 
]R+;L2(]R2)) < Be ; 

-;L4(r2)) < Ke 


Proof. The main idea is to use Lemma 2.1 with p = 4. As is divergence 
free and uniformly bounded, there exists K G Wq'‘^{A^) such that div/i^ = 
, so that 

divu^ = Vp^ 


-|- div — div K = 0. 


Using (2.5) in conjunction with (2.12), the estimates in Lemma 2.1 give: 


gl/2 

II^^IIl°°(L4) < CiSWVp^ ■ ^^11^00(^4) < QQ (l+^)/4 

de 


I E \\ 
\U 


which proves (0. Since uf belongs to 2))), we also obtain (0 

by Holder’s inequality: 


gl/2 
^1/2 




- ^4 (i+^)/ 4 C' 4 e||V())^||Lcx=||Ui ||L°°(L 4 (SuppV(ji=)). 
de 


We prove 0 in a similar way: 

eV2 


C 


||v/i44oo(i2) < C4-ty^j^c4v p^ ■ u 11^00(^4) < 

Oj£ de 


I E \\ 

\xi ||l°°• 
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Finally, to establish 0 , we observe that 

- u^Wl* < 11(1 - + Wh^Wi* < C 


M/2 




which concludes the proof. 


□ 


Remark 2.2. Since uf is hounded in L°°{L^) for any p < oo, Lemma 2.1 
with p = 2 also gives that 


%h"\\L^{L^) < C2£\\V<P^ ■ uf ||i^(i2) < 


U- 




where q = > 2 can be chosen as close as we want of 2. 

Nevertheless, for the sequel, we only need the case q = A (the case corre¬ 
sponding to the bound Q in Proposition 2-4 above). 


Remark 2.3. As we did not use (2.15) in the previous proof, we can give a 
more precise dependence of the constant Ke on the geometry and the Euler 
solutions (which are both fixed throughout). In fact, Ke depends only on 
K, a nd ||a;o||iinL°o. In particular, the bounds ([^, (Q in Proposition 
2.4 are linear in ||tt^||i:,oo. Hence, in these estimates Ke is of the form 


^\\^o\\HnL°°! where K depends only on 1C. 


3. Convergence of the initial velocity 


In this section, we discuss the convergence of the Navier-Stokes initial 
data (taken dependent on £ only) to the Euler initial data in the energy 
norm, as e goes to zero. 

Let wo G C“(M2 \ ([0,1] X {0})) if ^ G [0,1) and in ^“(M^ \ ([Q, 1]2)) if 
p = 1, then for s small enough, wq = 0 on for all i,j. For any e > 0 
there exists a unique vector held Uq, which is solution of (see e.g. [H] 
for a proof). We choose R so that Suppwo U C B{0,R). The function 

: z = xi ix 2 (iro)i(xi, X 2 ) - i(no) 2 (a;i, X 2 ) 

verihes the Cauchy-Riemann equations on B(0,Ry. Hence, is holomor- 
phic and admits a Laurent series decomposition fi{z) = YlT=i which 
allows us to conclude that Uq G n LF‘’°°{Bfi),RY) (where denotes 
the Marcinkiewicz weak L^-space). By elliptic regularity, Uq is bounded in 
B{0, R) n implying that 

(3.1) G nL2.~(Q£). 


This regularity of the initial data is sufficient to apply the result of Kozono 
and Yamazaki m. yielding existence and uniqueness of global solution to 
the Navier-Stokes equations (1.5). 

The goal of this section is to prove that Uq converges in uq dehned 

in (1.7). More precisely, we prove the following: 
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Proposition 3.1. Let QL he defined in (1.1)-(1.4) with d^ > e and let ujq G 
There exists a constant C, which depends only on K, such that, if 
lvq is supported in then 


11^0 ^o||l 2(Q£) < C||a;o||LinL“ 


e| lne| 
^(i+m)/2 ’ 


U,Q 


This proposition follows from the analysis developed in [6]. There, the 
authors have looked for the best condition on £ and d^ ensuring convergence 
of Uq to uq in LfiiOL), but unfortunately, they did not explicitly state the 
rate of convergence. To obtain the rate, we now briefly review the results in 
[6]. For brevity, we will denote Mq := ||wo||LinL°° ■ 


3.1. Correction and decomposition. We recall that uq has an explicit 
formula in terms of ojq via the Biot-Savart kernel (Equation ( |1.7[ )). However, 
there is no such formula available for Uq. In addition, uq is not tangent to 
and hence, it is necessary to introduce an explicit corrector, which is 
given in terms of: 

(1) a cutoff function equal to 1 close to ICfj: 

/x — zf,\ 

4>ij{x) :=</>( - - 1 ( with (j) given in 


( 2 . 2 )); 


(2) a biholomorphism T : M^\5(0,1) such that T( 2 ;) =/32: + /i( 2 ;) 

for some fi G M"*" and h a bounded holomorphic function. 

Then, the correction is defined by 


where 


:=^ [ In\x-y\u}o{y)dy 
2vr Jqs 




In- 


/3\x-y\ 








rMx)-TMyr\ 


In- 


(Xoiv) dy 

^o{y)dy, 


'n^ 


with 

/ X — zf 

Above, we have denoted by y* = the conjugate point to y across the 
unit circle in M^. This formula is related to the Biot-Savart law outside one 
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obstacle, and we can check that verifies the following properties: 

divn^ = 0 in • n = 0 on 


I 

JdK.j 


■ T ds = 0 for all i,j, lim |n^(x)| = 0. 


Then, we decompose uq — as 

(3.2) 

where 
wl{x) 


uo-v^ = J2wk, 

k=l 


I3\x - y\ 











In 


- V,jiy)\ 
\K,ix)-Tll^{yr\ 


'^o{y)dy, 


{x - y) 


|x - yp 


- {DTI,fix) 


^o{y) dy, 

{TMx)-VAy))- 


wlix) =—'^(j)lj{x){D7lA'^{x) 




ko(y) dy, 








nA\rfAA-%Ay)f 


^o{y)dy, 






where we have used that uiq is supported in 12^. In [6j and in [23] , the explicit 
formula for wf, is utilized to prove that converges to uq as e —>■ 0. In |6], a 
convergence analysis was also performed for <C e, which complicates the 
choice of the cutoff function. In our case where > s, the cutoff j defined 
above is sufficient to apply the estimates in |6]. We recall these estimates, 
adapted to the set-up of this paper, in the remainder of the section. 

3.2. Estimates on w, and w^. Following the proofs of Proposition 3.3 
and Proposition 3.4 in |6|, we first tackle estimates for wf and w^. When 
/C = 5(0,1), T = Id (so /3 = 1) and wf = = 0. In the general case, the 

main idea is to use that 


X — z: 




/3- 


X — z. 


= /3 Id 


VAf = ^ 

In |6l Proposition 3.3], it was proved that 


X — z, 






< CMqeI lne|, 


where C depends only on /C. Using that the supports of V(()f ^ are disjoint 
for distinct i,j (see (2.1)), we deduce that: 

< CMoe|lne|-||V(/>||Loc-^^^ = CMq 

de 
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In [H Proposition 3.4], it was also shown that, for all i,j and x G 11^, it 
holds 


f ( {x - y)^ 

L[\x-y\^ 




irf,ix)-Tf,iy)) 


'6 


< CMo 


wo(y) dy 


.1/2 


+ 


\x — z: 




where C depends only on fC. As there exists <5 > 0 such that i?(0,(I) C JC, 
Supp(^|^- C B{zlj,2y/2£) \ B{zfj,d£), and we can estimate the norm of 
as follows: 




e 

^(i+m)/2 


+ CMo£ 




' B(0,2V2e)\B{0,Se) 


dx \l/2 


< CMo 


e| lne|^/^ 
^(i+/^)/2 • 


3.3. Estimates on tcl and tcf. When K, = i?(0, 1), we have T^j{y)* = 

y — ■ 

We observe that, even though \y — zf A can be of order e (this 

\y-zlj\^ 

is, in fact, the case if y G dB{zf j,e)), Tfj{y)* is still small compare to 7lj{x), 
since \Tfj{x)\ > 1. Hence, rcl and wl should be small. 

More precisely, in [HI Proposition 3.5], it was proved that, for all i,j, 


< CMo£ 

L°°(Supp V</? j.) 

where C depends only on 1C. Therefore, we find that 

lhillL2(0-) < CMo£-\\V^L-- 

Cte Cte 

Concerning the estimate on u;|, the proof of Proposition 3.6 in [6] shows 
that, for all i,j and x G 


In 


\T-\ 

' * j 


X 


- T- 

1,3 


{yr 




\T^- 


(x)| 


-wo(y) dy 



vA^)-‘n:jiyr 

%jix) - Tfj{y)*\^ 



wo(y) dy 


S + 




where C depends only on JC. Therefore, after the change of variable = 
eBfjix), we can estimate the norm in the same way as for (see [6] for 
details), and we conclude that 


||i«4|Il2(Oe) < CMo 


e| lne|^/^ 
^(i+m)/2 ■ 
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3.4. End of the proof of Proposition 3.1| . Collecting the estimates on 
i(;| for A: = 1,..., 4 gives that for any e > 0, there exists a constant C > 0 
such that 

Cie 

Finally, we observe that Uq — is the Leray projection of uq — (see |6l 
Sect. 3.2]). Hence, the L^-orthogonality of this projection implies 

- < ho - ^^*"1^2(0-) < 

dp 


We conclude the proof of Proposition 3.1 by the triangle inequality. 


4. Energy estimate 

In this section, we establish convergence of the solution u'^’^ of the Navier- 
Stokes equations in to the solution of the Euler equations in the whole 
plane, provided H u for some given positive constant A. The proof 

is performed by an energy estimate on the difference — u^, where 


2.4). This strategy is 


vanishes on and is close to (see Proposition 
related to the proof of Kato’s criterion (see HZlISlEe]) and was used in m 
to treat the case of one shrinking obstacle. 

Heuristically, it is possible to pass to the zero-viscosity limit even in the 
presence of the boundary layer due to the no-slip boundary condition, since 
under the conditions of the theorem below the contribution to the energy 
from the Euler solution restricted to the obstacles is negligible. Nevertheless, 
in the case of infinitely-many obstacles a dependence on dg, which reduces 


the expected rate of convergence of ^/u, valid for one shrinking obstacle, ap¬ 


pears in the convergence estimate (4.2). This dependence can be interpreted 
as a “ghost” of the perforated domain. 


Theorem 4.1. Let he the solution of the Euler equations (1.6) on 
[0,oo) X with initial condition uq (given in (1.7) in terms of an initial 
vorticity ujq G C'“(M^)j. Let he the solution of the Navier-Stokes prob¬ 


lem (1.5) on [0, oo) X (where is defined in (1.1)-(1.4) with d^ > e) with 
initial velocity Uq^ G L°°r\L'^’°°{Ll^), satisfying the divergence-free condition 
and no-penetration boundary condition, and such that — uq G LP‘{LI‘^). 
Then, there exists a constant A, depending only on JC, such that if 

Av 


(4.1) 

then for any 0 < T < oo, 


< 


^(i+n)/2 IjwollLinLo 


(4.2) sup - u^\\L 2 t^e) < Bt 

0<t<T 


^( 1 + m )/2 + IIV -^ 0 || l 2 ( O ^) 


where Bt is a constant depending only on T, ||a;o||Liniunoo, and 1C. 
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We recall that depends on the conhguration of the obstacles, which 
are arranged according to the parameter fi G [0,1]. 

Remark 4.1. Before giving the proof of the theorem, we list two immediate 
corollaries: 

(a) By choosing as initial data for Navier-Stokes Uq^ = Uq, where Uq 
solves (1.9), Proposition 3.1 and Theorem 4-1 immediately give that 


sup \\u 
0<t<T 


Eli D ( e|lne| \ 

- U 11^2(0.) < Bt + ^( l +/.)/2 j ’ 


where 

e\ lne| < Cy/e < C^/u. 

This estimate then implies Theorem 
(b) In the case of an arbitrary, but fixed and independent of e, num¬ 
ber of shrinking obstacles, we recover the rate of convergence of y/v 
obtained for one obstacle in [T6] . 


Proof of Theorem \4.1\ With defined in Proposition 2.4 we introduce 

:= 

which is divergence free and verifies the Dirichlet boundary condition on 
We observe that satisfies: 

uf = (jf uf — hi = fp{—u^V— Vp®) — /if. 

Hence, satisfies: 

W'f'^ - = - Vp’'’^ + uAu^ - 

+ ■V)u^ + Vp^] + hl. 

Multiplying this equation by and integrating by parts, we obtain: 


2dt" 


= — u 


: Vu^ dx- dx 


(4.3) 


+ / ■ V)n®] dx + / 




+ 


W‘'’^ • /if dx. 




Under the assumption on the initial data, u'^’^ and have enough regularity 
to justify the integration by parts and, hence, the energy identity. 

We group the trilinear terms in (4.3) together and tackle each of the other 
terms separately. We start with the trilinear terms: 


X := - 




• [((ii'^’" • V)ii^’") - • V)u^)] dx, 
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which we rewrite, owing to and the divergence-free condition 

on as 


X= - 




^ 3 . _ / ^ru,e . 




+ / • V)u^] dx. 

Jn^ 

We will ntilize that is close to both (j)‘" and to control the last two 
terms. Hence, we add and subtract • [(n^ • V)n^] dx: 

I=- [ • V)n^] dx - [ • [{u^ ■ V){u^ - n®)] dx 




+ / ■ V)u^] dx 


in^ 


= - [ • [(W^’= • V)n^] dx+ f (n= - n®) • [(n= • V)W''’^] dx 


(4.4) 


+ / W^'’^ ■ [{h^ ■ V)u^] dx, 


where we have integrated by parts in the second term and have used that 
4>^u^ — . 

We bound each of the three integrals in ( |4.4[ ) above separately. We set: 
Ii :=- [ • V)n^] dx. 

We have 

Vrt^ = (P^Vu^ + V(/>^ ® - V/i^ 

with the last two terms supported on (defined in (|2.4[)). Hence, using esti¬ 


mates (2.5) on Vcj)^^, (2.12)-(2.15) on , and estimate (l3|) in Proposition 
on V/i^, we have thanks to Lemma |2.3[ 

+ ll^"’^lli4(A.)llvh^llL2 

(4.5) 


2.4 


where depends only on uq and /C. From the discussion in Remark 2.3 


it's can be taken of the form K 3 = K^\\ojQ\\iif^ioo, where depends only 
on /C. For 

X 2 := [ {u^ - u^) ■ [{u^ ■ V)W’^] dx, 

Jn^ 

we note that supp(u^ —n®) C Ag. Then, Holder’s and Cauchy’s inequalities, 
together with estimates ( |2.5[ ) on 1 — cp^, (2.12) on , and estimates 0 and 
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Q in Proposition 2.4, give that 






d. 


(4.6) + 

where K 4 depends only on wq and 1C. 

Lastly, we bound 

X3 := f . [{h^ . V)u-®] dx. 

Using that supp C we apply Holder’s and Cauchy’s inequalities 
again, together with estimates Q of Proposition 2.4 for h^, (2.13) for Vu^, 
and Lemma 2.3 to obtain: 

IX3I < ||lU"’^||i4(^^)||h^||i4||Vn®|U2 

(is 

(4.7) + 

where depends only on ojq and 1C. 

From (^) and (|4.5), ( |4.6[ ), ( |4.7[ ), it follows that 

(4.8) 


|2:| < (4 




+ Ke 




where iLg depends only on wq and 1C. 

We now turn to the remaining integral terms in (4.3). We begin with 

J ■=-,,[ :Vu^dx. 

Applying Cauchy-Schw artz followed by Cauchy’s inequality, owing to esti¬ 
mate ([^ in Proposition 2.4 for V/i^, (2.5) for and (2.12)-(2.13) for 
gives: 


\J\ < v^l|vw"’^||i2(o.) 

Eu ^*2 


(4.9) 

where K^ depends only on ojq and JC. 


Ke y 
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Finally, we consider the last two integral terms in (4.3). We set: 

■= [ • Vp® dx = - [ ■ V<i)^ dx. 






Consequently, by Holder’s and Cauchy’s inequalities, and the Poincare’s 
inequality (2.11), using that V(/>^ is supported on together with estimates 
(2.5) for V(f)^ and ( 2.14| ) for p®, we have 

I'Hil < \\W''’^\\L2i^A,)\\'^^^\\L^{A^)\\P^\\L^iA,) 

C2 


ds 


(4.10) 




where depends only on ojq and 1C. 
We tackle the last integral term 


^2-.= - W^'^-hldx. 
Jn^ 


As /if is supported on A^, a similar application of the Cauchy-Schwartz, 
Cauchy’s inequalities, and the Poincare’s inequality (2.11), together with 
estimate ([^ in Proposition 2.4 on /if, gives: 


(4.11) 


de 

.-3 

7-u,eii2 




where Kg depends only on cog and 1C. 


Putting together the bounds (4.8) obtained for X, (4.9) for J", and (4.10)- 


(|4.11|) for F/, yields the following a priori energy estimate: 


2dt" 


3u e 


/3u 

-VT”^ ^^(1+m)/2 




+ Coe^ 


C'oq|^i^,£||2_ 


lL2(ne) 


+ 


'dl+^ 


Kg + Kg + Kg e" 


+ Kj 




We will apply Gronwall’s Lemma to the previous inequality. We therefore 
need the coefficient of ||V 1 F ^’^||^2 to be strictly less than u, which forces the 
following relation between u, e and d^: 

^( 1 +/ 4)/2 - 

where A = 1/(4A'3). Recalling that K^ = K^WuJoWiif^ioa, we obtain the 
compatibility condition (4.1) with A = l/(4iF3). This condition also implies 
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that 




(^6 + Ks + Kq e34i+M)/2^ < Kwiy, 


where Kiq depends only on cjq and JC. 

For any given 0 < T < oo, we next define 

Kt := Coe^o'^, 


with Co as in (2.15). Then, by Gronwall’s Lemma we have that, for any 
t G [0, T], 




Kiqv Kji' 


+ 


Kt Krd, 


i+/i 


o'^Krt 


which implies 
(4.13) 


Ilir^-'lo.oilLw + y^ 


Kiqv 

iLj' 


+ 


I Kri^ 

Krdl^'^ 


^Krt 


To conclude the proof, we apply the triangle inequality to obtain 

- n®)(t, 0^2(0^) < ||VF"’"(L + ll(«" - 

and 

\\W'^’^{0, •)||l2(qe) < ||Uq'^ — tto||L2(Qs) + ||u^(0, •) — Mo||L2(f2£)- 

It then follows from (2.5), (2.12), and point 0 in Proposition |2.4| that 
||(«^ - U^){t, •)llL2(n=) < ||?x®||l<-|| 1 - + ||/i^lU4||l||i4(^^) 


(4.14) 


< c 


41 +M )/2 


< cv, 


for all t G [0,oo), using also condition (4.12). Lastly, we combine (4.13)- 
(4.14) and note that 

fv 


I/, yjv < C 


41 +m )/2 ’ 


for some C > 0, as Z 2 and dg are both bounded above. Estimate (4.2) then 
follows. 

□ 


Remark 4 . 2 . The dependence on d^ in estimate (4.2) comes from the con¬ 
tribution of J, more precisely from the terms: 

V f : [V(/>^ (g) li®] dx and v f : VK dx, 

Ja^ Ja^ 

Without these terms, we would establish a rate of convergence of order y/v 
instead to in the case of a fixed number of inclusions. Even 

though these integrals are taken on A^, the measure of which is small, the 
L2 -norm of is unbounded ( |2.5 ), and hence we are not able to prove 
these terms are of the same order as the other terms in the energy estimate. 
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To improve the energy estimate, we can ask whether there exists a better 
choice for (ff than our construction (2.2)-(2.3). Namely, we look for a 
that vanishes in B{{),deY, equal to 1 in B{0,e), and such that ||V ^^||£,2 
is the smallest possible. This optimization problem has a unique solution, 
which is the solution of 

Acj)^ = 0, ^^\dB{0,e) = 1) <^^laR(0,4) = 0- 

This slution has an explicit expression of the form: 


= 


Mg) 

Mf)' 


Defining 
(4.15) 
it follows for de s. 


(p%x) {x- zN) , 

i=i j=i 


||V# 


Il2 ^ 


< c. 




< 


c 


Inel 


which is a slight improvement over (|2.5|). We cannot expect a rate better than 


( l +^)/2 


4 

term’ 


-v/jlnej 


in (4.2), because, at v fixed, it is well known that a “strange 


appears when 

d(Y+>^)D^\lne\ ^ C > 0. 

(See [U [8] for the case p. = I and [2] for the case p = 0.) 

We note that for p = 0, the situation in which the distance is smaller than 
the critical distance (namely, y^l lne| —)• Oj was not treated in [2] at fixed 
viscosity. This case was considered only when the inclusions are distributed 
in both directions (related to the case p = 1) filling a bounded domain [T]. 
Indeed, this set up implies a uniform Poincare estimate. Allaire investigated 
the limit of /{d^y/l lne|), which in turn implies a limit for the term: 


\Vu''Y\'^(f^\dx = 


f 



dsV\ 




dx. 


Therefore, guided by the optimal results for the Laplace and Stokes equa¬ 
tions, it would be natural to consider cj)'^ defined as in (4.15) rather than in 


(2.3). However, the more natural construction of the cut-off function leads 


to a series of technical complications. Besides the adaptation of (2.5), we 


have used several times elliptic estimates on a domain of the form U = 
(—2, 2)^ \/C and the fact that + sU). On the other hand, dilating 

B{t),de) \ B{0,e) by a factor of e results in a set the size of which grows 
unboundedly os e —)> 0. To adapt the approach of Section\^ to the improved 
setting, we would therefore have to: 
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estimate the Bogovskii constant Cp in Lemma 


2.1 


in 


terms of e, 


which should be uniformly bounded thanks to [9l Theorem III.3.1]; 
estimate the Poincare constant Ki in Lemma\2.1^ 


estimate the embedding constant K in the proof of Lemma 2.3. 
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